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Abstract 



By exploiting the mathematical analogy between the propagation of sound 

in a non-homogeneous potential flow and the propagation of a scalar field in a 

background gravitational field, various wave "energy" and wave "momentum" 

conservation laws are established in a systematic manner. In particular the 

acoustic energy conservation law due to Blokhintsev appears as the result of 

the conservation of a mixed co- and contravariant energy- momentum tensor, 

while the exchange of relative energy between the wave and the mean flow 

mediated by the radiation stress tensor, first noted by Longuet-Higgins and 

Stewart in the context of ocean waves, appears as the covariant conservation 

of the doubly contravariant form of the same energy- momentum tensor. 
PACS numbers: 43.28.Py, 43.20.Wd, 43.25.Qp, 67.40.Mj 



1 



Typeset using REVT[t;X 



I. INTRODUCTION 



Many discussions of the "energy" and "momentum" associated with waves propagating 
through moving fluids can be found in the physics engineering 0-^, and mathematical 
fluid mechanics literature |3"16|- Various definitions are proposed, some of which lead to 
conserved quantities, and some to quantities that are not conserved but instead exchanged 
between the wave and the mean flow. In part the multiplicity of definitions is due to difficulty 
in deciding what fraction of the energy or momentum of the system properly belongs to the 
wave and what fraction should be associated with the moving medium. It is also often 
unclear how to divide equations expressing conservation laws into terms relating to the 
conserved quantity, and terms acting as sources for this quantity. Related to these primarily 
cosmetic problems are more fundamental issues as to whether the "energy" or "momentum" 
under discussion is the true newtonian energy or momentum, or instead pseudo-energy and 
pseudomomentum. Thus we have the question "what is the momentum of a sound wave" 
raised by Sir Rudolf Peierls in his book Surprises in Theoretical Physics fl^, and the salutary 
polemic "On the 'Wave Momentum' Myth" by M. E. Mclntyre 



The most extensive analyses of conserved wave properties have been carried out by the 
fluid mechanics community p|-p!5|. Typically these papers adopt a lagrangian (following 
individual particles in the flow) or mixed lagrangian-eulerian approach, as opposed to the 
purely eulerian (describing the flow in terms of a velocity field) approach which would be 
most familiar to a physicist. In addition a physicist reading this literature feels the lack of 
a general organising principle behind the definition and derivation of the conservation laws. 
The present paper is intended to remedy some of these problems — at least for the special 
case of sound waves propagating through an irrotational homentropic flow. Although rather 
a restricted class of motions, this is still one of considerable interest in condensed matter 
physics. It includes phonon propagation in a bose condensate, and therefore lies at the heart 
of the two fluid model of superfluids. By exploiting W. Unruh's ingenious identification 
I9| , |20[] of the wave equation for sound waves in such a flow with the equation for a scalar 
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field propagating in a background gravitational field, I extract the conservation laws from 
the principle of general covariance. Deriving the conservation laws in this way may seem 
like a case of taking a sledge-hammer to crack a nut, but the formalism is familar to most 
physicists, automatic in application, and the ambiguities in defining the conserved quantities 
turn out to lie in the choice of whether to identify the energy-momentum tensor as T^*^ or as 
T'^^. Also, when quantities are not conserved, as is the case of the wave momentum in a shear 
flow, their sources arise naturally from the connection terms in the covariant derivative. 

In section two I discuss the action describing the irrotational motion of an homentropic 
fluid. In section three I derive Unruh's equation from the action principle. In section four I 
explain why we often need information beyond the solutions of linearized wave equation, and 
in section five derive the conservation equations that follow from the linearized equation. 
Section six interprets these equations in terms of the motion of phonons. The discussion 



section repeats the warnings from [|T^ that the appearance of quantum quasi-particles in any 
argument is a sure sign that the we are considering pseudomomentum and not the actual 
momentum of the system. 

The work reported here was motivated by a desire to understand better the role of 
acoustic radiation stress in the two-fluid model of a superfluid. It may be relevent to the 



recent controversy p2|- P^] over the lordanskii force acting on a vortex moving with respect to 
the normal fluid component. The use of the Unruh formalism in this context was suggested 



by Volovik 25 



II. THE ACTION PRICIPLE 

The most straight forward way of deriving conservation laws starts from an action prin- 
ciple. Noether's theorem then provides us with an explicit formula for a conserved quantity 
corresponding to each symmetry of the action. In fluid mechanics, unfortunately, at least 
when we restrict ourselves to an eulerian decription of the flow field, action principles are 
in short supply. Of course there must exist some action principle because ultimately the 
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fluid can be treated as a system of particles. A particle-based action, however, requires a 
lagrangian description of the flow. When it is re-expressed in eulerian terms constraints 
appear, and these limit its utility. 

If we restrict ourselves to flows that are both irrotational and homentropic — the latter 
term meaning in practice that we assume the pressure to be a function of the fluid density 
only — then the number of degrees of freedom available to the fluid is dramatically reduced. 
In this case the eulerian equations of motion are derivable from the action ||2^ 



S = jdS {p0 + \p{Vct>f + n(p)} . (2.1) 

Here p is the mass density, (f) the velocity potential, and the overdot denotes differentiation 
with respect to time. The function u may be identified with the internal energy density. 
Equating to zero the variation of S with respect to yields the continuity equation 

P + V-(H = 0, (2.2) 

where v = V0. Varying p gives Bernoulli's equation 

0+ V + /i(p) = 0, (2.3) 

where p{p) = du/dp. In most applications p would be identified with the specific enthalpy. 
For a superfluid condensate the entropy density, s, is identically zero and p is the local 
chemical potential. 

It is worth noting that our action could not have arisen from some rewriting of the 
action for the motion of a system of individual particles. We are allowing variations of p 
without requiring simultaneous variations of 0, and such variations conjure new matter out 
of nothing. 

The gradient of the Bernoulli equation is Euler's equation of motion for the fluid. Com- 
bining this with the continuity equation yields a momentum conservation law 

dt{pvi) + dj{pVjVi) + pdip = 0. (2.4) 
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We simplify ( |2.4| ) by introducing the pressure, P, which is related to /i by -P(p) = / prf/U. 
Then we can write 

dt{pvi) + d,Ii,i = (2.5) 

where Iljj is given by 

liij = pViVj + 6ijP. (2.6) 

This is the usual form of the momentum flux tensor in fluid mechanics. 

The relations fi = du/dp and p = dP/dp show that P and u are related by a Legendre 
transformation, P = pp — u{p). From this and the Bernoulli equation we see that the 
pressure is equal to minus the action density, 

-P = p0+^p(V0)2 + u(p). (2.7) 

Consequently we can write 

n,, = pa,0a,-0 - 5,, |p0 + ^p(V0)2 + m(p)| . (2.8) 

This is the flux tensor that would appear were we to use Noether's theorem to derive a law 
of momentum conservation directly from the invariance of the action under the translation 
0(r) 0(r — a), p(r) p{r — a.). This is not a trivial point because there are at two similar, 
but distinct, notions of "momentum". True momentum is associated with the symmetry of 
the action under a simultaneous translation all the particles in the system. Its conservation 
requires an absence of external forces. Pseudomomentum is the quantity that is pre- 
served when the action is left invariant when the disturbance in the medium is relocated, 
but the reference position of each individual particle is left unchanged. Conservation of 
pseudomomentum requires homogeneity of the medium rather than of space. Replacing the 
field 0(r) by 0(r — a) would normally correspond to the latter symmetry, but, because of 
the absence of explicit particles, at this point in our discussion the two concepts coincide. 
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III. THE UNRUH METRIC 



We now obtain the linearized wave equation for the propagation of sound waves in a 
background mean flow. Let 

(p = 0(0) + 0(1) 

P = P(o)+P(i)- (3.1) 

Here 0(o) and p(o) define the mean flow. We assume that they obey the equations of mo- 
tion. The quantities 0(i) and p(i) represent small amplitude perturbations. Expanding S to 
quadratic order in these perturbations gives 

S = So + J c^'x|p(i)0(i) + ^-^^p?i) + ^P(o)(V0(i))2 + P(i)V. V0(i)|. (3.2) 

Here v = V(o) = V0(o). The speed of sound, c, is defined by 



P(o) dp 



(3.3) 



P{0) 



or more familiarly 



o dP , , 



The terms linear in the perturbations vanish because of our assumption that the zeroth-order 
variables obey the equation of motion. 

The equation of motion for p(i) derived from ( p.2| ) is 

In general we are not allowed to substitute a consequence of an equation of motion back 
into the action integral. Here however, because p(i) occurs quadratically, we may use ( p.5|) 
to eliminate it and obtain an effective action for the potential 0(i) only 

^(2) = / d'x {ip(o)(V0(i))2 - ^(0(1) + V ■ V0(i))^} . (3.6) 



The resultant equation of motion for 0(i) is [IS, 2 



(I + V ■ v) ^ (I + V . v) 0(1) = V(p(o)V0(,)). (3.7) 

Note that in deriving this equation we have not assumed that the background flow v is 
steady, only that it satisfies the equations of motion. Naturally, in order for our waves to 
be distinguishable from the background flow, the latter should be slowly changing and have 
a longer length scale than the wave motion. 

Now ( |3.7D is equivalent to a perhaps more familiar equation 

(I + ^ ■ ^) ? (s + ' ■ ^) = 

as can be seen by using the mass conservation equation dtP(Q) + V ■ P(o)V = 0, but the form 
( p.7|) has the advantage that it can be written asQ 

^5/.v^^7^"5.</'(i) = 0, (3.9) 

where 

^9r='-f('' (3.10) 
^ \ V, vv^ — c 1 / 

This is perhaps most easily seen by observing that the action ( ^.6| ) is equal to — S" where 

S = j d^x iv^^'^'^a^0(i)9,</>(i) = j d'^x v^L. (3.11) 

Equation (|3.9|) has the same form as that of a scalar wave propagating in a gravitational 
field with Riemann metric The idea of writing the acoustic wave equation in this way. 



as well as the general relativity analogy, is due to Unruh |1T9| , |20|1 . I will therefore refer to (7^1/ 
as the Unruh metric. 

The 4-volume measure is equal to pfo^/c, and the covariant components of the 

metric are 



use the convention that greek letters run over four space-time indices 0, 1, 2, 3 with = t, while 
roman indices refer to the three space components. 
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9,^ = ^-fy ^ .]■ (3-12) 

The associated space-time interval is therefore 

ds^ = ^ [c^dt^ - Sij{dx' - v'dt){dx^ - vHt)] . (3.13) 

Metrics of this form, although without the overall conformal factor p{o)/c, appear in the 
Arnowitt-Deser-Misner (ADM) formalism of general relativity |^. There c and — are 
refered to as the lapse function and shift vector repectively. They serve to glue successive 
three-dimensional time slices together to form a four dimensional space-time I^H]. In our 
present case, provided P(o)/c can be regarded as a constant, each 3-space is ordinary flat 
equipped with the rectangular cartesian metric Qij^^^^^ = 6ij — but the resultant space-time 
is in general curved, the curvature depending on the degree of inhomogeneity of the mean 
flow V. 

In the geometric acoustics limit sound will travel along the null geodesies defined by g^^- 
Even in the presence of spatially varying p(o) we would expect the ray paths to depend only 
on the local values of c and v, so it is perhaps a bit surprising to see the density entering 
the expression for the Unruh metric. However an overall conformal factor does not affect 
null geodesies, and thus variations in p(o) do not influence the ray tracing. For steady flow, 
and in the case that only v is varying, it is shown in the appendix that the null geodesies 
coincide with the ray paths obtained by applying Hamilton's equations for rays 

to the appropriate Doppler shifted frequency 

w(x,k) = c|k| +v-k. (3.15) 

When V is in the x direction only, we can also rewrite ds"^ as 

ds'^ = M |_ (^dx -{v + c)dt) {dx -{v- c)dt) - dy^ - dz^] . (3.16) 

This shows that the x — t plane null geodesies coincide with the expected characteristics of 
the wave equation in the background flow. 
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IV. MOMENTUM FLUX 



The fluid in a sound wave has average velocity zero, but since the fluid is compressed 
in the half cycle when it is moving in the direction of propagation and rarefied when it is 
moving backwards there is a net mass current (and hence a momentum density) which is 
of second order in the sound wave amplitude Qq. This becomes clearer if one solves the 
equation 



for the trajectory x{t) of a fluid particle. This equation is non- linear (x appears inside the 
cosine), and solving perturbatively one finds a secular drift at second order in qq. 



Although the time average of the eulerian fiuid velocity, v, is zero, the time average of the 
lagrangian velocity, vl — x, is not. The difference beweeen the two average velocities is the 
Stokes drift. The Stokes drift is 0{al) while the wave equation is accurate only to 0{ao), so 
care is necessary before using its solutions to evaluate the mass current. Similar problems 
occur in defining the energy density and energy and momentum fiuxes which also require 
second order accuracy. 

We can expand the velocity field as 



where the second-order correction V(2) arises as as consequence of the nonlinearities in the 
equations of motion. This correction will possess both oscillating and steady components. 
The oscillatory components arise because a strictly harmonic wave with frequency ujq will 
gradually develop higher frequency components due to the progressive distortion of the wave 
as it propagates. (A plane wave eventually degenerates into a sequence of shocks.). These 
distortions are usually not significant in considerations of energy and momentum balance. 




(4.1) 




(4.2) 



V = V + V(i) + V(2) H , 



(4.3) 
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The steady terms, however, represent O(ao) alterations to the mean flow caused by the sound 
waves, and these often possess energy and momentum comparable to that of the sound field. 

Even if we temporarily ignore these effects and retain only V(i) as determined from the 
linearized wave equation, the density and pressure will still have expansions 

P = P(o) + P(i) + P(2) H 

P = P(o) + P(i) + P(2) + • • • . (4.4) 

As before, the grading (n) refers to the number of powers of the sound wave amplitude in an 
expression. The small parameter in these expansions is the Mach number given by a typical 
value of i'(i) divided by the local speed of sound. 

Consider for example the momentum density pv and the momentum flux 

Tlij = pViVj + SijP. (4.5) 

It is reasonable to define the momentum density and the momentum flux tensor associated 
with the sound field as the second order terms 

j^'''°'^°'^^ = (pv) = (p(i)V(i)) + v(p(2)), (4.6) 

and 

j^(phonon) ^ ^^^^ {v^i)iV^i)j) + Vi{p^i)V^i)j) + Vj {p^i)V^i)i) + 6ij{P^2)) + (P(2) ) • (4.7) 

(The angular brackets indicate that we should take a time average over a sound wave period. 
There is no need to consider terms first order in the amplitude because these average to zero.) 
We see that to compute them we need to consider the second order contributions to both P 
and p. 

We can compute P(2) in terms of first order quantities from 
and Bernoulli's equation in the form 
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A/i = -0(1) - ^(V0(i))' - V ■ V0(i), (4.9) 

together with 

^ = p fZ = ^=P_ (4 10) 

Expanding out and grouping terms of appropriate orders gives 

^(1) = -P(O)(0(1) + V ■ V0(i)) = cV(i), (4.11) 

which we already knew, and 

P(2) = -P(o)^(V</>(i))2 + ^^(0(1) + V • V</.(i))^ (4.12) 

We see that P(2) = y/—gL where L is the Lagrangian density for our sound wave equation. 

To extract p(2) in this manner we need more information about the equation of state 
of the fluid than is used in the linearized theory. This information is most conveniently 
parameterized by the logarithmic derivative of the speed of sound with pressure (a fluid- 
state physics analogue of the Griineisen parameter). Using this together with the previous 
results for P(2), we find that 



1 „ 1 2 d,\nc 



p(2) = -^Pm - — P(i) 



P(o) ^ ' d\\ip 



(4.13) 

P(0) 



V. CONSERVATION LAWS 

While we cannot compute the "true" energy and momentum densities and fluxes without 
including non-linear corrections to the motion, it is often more useful find closely related 
quantities whose conservation laws are a consequence of the linearized wave equation, and 
which therefore provide information about the solutions of this equation. Our "general 
relativistic" formalism provides a sytematic way of finding such conserved quantities. It is 
well known PS| that any action S automatically provides us with a covariantly conserved 



and symmetric energy-momentum tensor T^^, defined by 
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The functional derivative is here defined by 



SS = Jd'xV^j^6g^^. (5.2) 

It follows from the equations of motion derived from 5* that 

D^T^" = 0, (5.3) 

where is the covariant derivative. For example 

D^A'^\ = d^A^\ + r^^^A''^ + Tl^A^\ - Tl^A^'^^. (5.4) 

The are the components of the Levi-Civita connection compatable with the Unruh 
metric, viz. 

r^, = ^7"1/57,p], (5.5) 

where 



For our scalar field 

rpflU 



S^^djS'^^d) - g'^'^ . (5.7) 



The derivatives with raised indices in ( p. 71 ) are defined by 



9V(i) = g^-d.cP^i-) = ^(0(1) + V ■ V0(i)), (5.8) 

P(0)C 



and 



= g'^d^^ii) = (t;.(0(i) + V ■ V0(i)) - 0^9,0(1)) . (5.9) 

P(0)C ^ ^ 

Thus 
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P(0) V C2 

= ^Pm- (5.10) 

P(0) 

The last two equalities serve as a definition of Wr and p(2). The quantity PFr is often 
decribed as the acoustic energy density relative to the frame moving with the local fiuid 
velocity [|r^. Because its conservation law will depend on the steadiness of the fiow rather 



than the absence of time-dependent external forces, it is more correctly a pseudo-energy 
density. 

Using ([4.11|) , and ( [4.12|) in the form 



= ^P(2), (5.11) 
^ P(0) 

we can express the other components of ( p. 71 ) in terms of physical quantities. We find that 



P(0) 

= ^(P(i)%)* + ^*P(2)j- (5-12) 

The first line in this expression shows that, up to an overall factor, is the energy flux - 
the first term being the rate of working by a fiuid element on its neigbour, and the second 
the advected energy. The second line is written so as to suggest the usual relativistic iden- 
tification of (energy-fiux) /c^ with the density of momentum. This interpretation, however, 
requires that p(2) be the second order correction to the density, which it is not. 
Similarly 

T'^ = — (p{o)^(i)i^^(i)j + ^^iP(i)^(i)i + VjP{i)V{i)i + SijP(2) + ViVjP(2)) ■ (5.13) 

P(0) 

We again see that if we identify p(2) with p{2) then T'-' has the exactly the form as we expect 
for the second order momentum fiux tensor. 

The reason why the linear theory makes the erroneous identification of p(2) with p(2) is 
best seen if we set v = const. Then the equation 
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dtT°^ + d^T^ = 0, (5.14) 

holds. This reads 

{dtp{2) + di{p{i)V(i)i + ViP(2))) = 0, (5.15) 
and looks very much like the second order continuity equation 

dtp{2) + di{v{2)p{Q) + P(i)^(i)i + Vip{2)) = 0, (5.16) 

once we ignore V(2). When we go beyond the linear theory ( |5.1(j| ) provides a source or sink 
term in the mass conservation equation for (v(2)) and is not an equation determining 

P(2)- 

We can also write the mixed co- and contra- variant components of the energy momentum 
tensor T^^^ = T^^gx^ in terms of physical quantities. This mixed tensor turns out to be more 
useful than the doubly contravariant tensor. Because we no longer enforce a symmetry 
between the indices p and u, the quantity Wr is no longer required to perform double duty 
as both an energy and a density. We find 

T'o = -^{Wr+ P(1)V(1) ■ v) 

^0 = ^ — + V ■ V(i) ) (p(o)t^(i)i + P(i)V(o)i), (5.17) 
P(o) \P(o) / 

and 

P{0) 

T) = --^ {p{o)V(i)iVii)j + t^iP(i)t^(i)i + 5ij^(2)) . (5.18) 
P(o) 

We see that p(2) does not appear here, and all these terms may be identified with physical 
quantities which are reliably computed from solutions of the linearized wave equation. 

The covariant conservation law can be written as either D^T^'^ = or as D^T^^ = 0. The 
two equations are consistent with each other because the covariant derivative is defined so 
thatZ^Afl'^ti^ = g^lvD\. To extract the physical meaning of these equations we need to evaluate 
the the connection forms V^^. 

14 



In what follows I will consider only a steady background flow, and further one for which 
Po, c, and hence = P(o)/c can be treated as constant. To increase the readabilty of 

some expressions I will also choose units so that po and c become unity and no longer appear 
as overall factors in the metric or the four dimensional energy-momentum tensors. I will 
however reintroduce them when they are required for dimensional correctness in expressions 
such as P(o)V(i) or Wr/c^. 

From the Unruh metric we find 

[tj,k]=0 

[^J.O] = -(diVj + djVi) 
[^0, j] = ^{diVj - djVi) 
[0^,0] = [zO,0] = -^9i|t;p 
[00, z] = ]^di\v\^ 

[00,0]=0. (5.19) 

I have retained the expression \{diVj — djVi) in [iO,j], since it is possible that our wave 
equation has greater generality than its derivation. 
We therefore find 

rSo = ^(v-V)|t;r 

1 1 

r°o = -i^di \v\^ + -Vj{d,vj - djVi) 



■'I I 

^)k = \<d3'^k^dkVi). (5.20) 



2 

From ( [5.6| ) we have 

1 d 



-g dx^ 

so, with ^/—g = const., the trace r|^^ is zero. One may verify that the above expressions for 
V^^ obey this identity. 
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We now evaluate 




(5.22) 



After a little algebra we find 



(5.23) 



Note the non-appearence of p(i) and p(2) in the final expression — even though both quan- 
tities appear in T^'^. 

The conservation law therefore becomes 



This is an example of the general form of energy law derived by Longuet-Higgins and Stuart, 
originally in the context of ocean waves |T^,|T3[. (See also |^ for a slightly earlier, but less 
general, case.) The relative energy density, Wr = T^^, is not conserved. Instead an observer 
moving with the fluid sees the waves acquiring energy from the mean flow at a rate given by 
the product of a radiation stress Sjj with the mean-flow rate of strain. Such non-conservation 
is not surprising. Seen from the viewpoint of the moving frame the flow is no longer steady, 
while (pseudo) energy conservation requires a time-independent medium. 

Notice that, since we are assuming that p(o) is a constant, we should for consistency 
require V ■ v = 0. Thus the isotropic part of the radiation stress (the part oc 6ij) does 
no work. This is fortunate because the non-linear theory shows that the isotropic radiation 
stress contains a part dependent on dlnc/dln p which is missed by the linear approximation, 
(see however, pT| ) 

We now examine the energy conservation law coming from the zeroth component of the 
mixed energy-momentum tensor. We need 




(5.24) 



where 



(5.25) 
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= d;r% - [zO, 0]T'° - [00, i]T^' - [iO, (5.26) 

We now observe that = while [00, i] = -[iO,0], and that [iO,j] = -[jO,i], while 
rpij _ rpji^ Thus the connection contribution vanishes. This form of the energy conservation 
law is therefore 

dt (Wr + P(i)V(i) -v) +di I (— + V ■ V(i))(p(o)t;(i)i + P(i)t;(o)i) ) = 0- (5.27) 
^ ^ V P(o) / 

Here we see that the combination Wr + P(i)V(i) ■ v does correspond to a conserved energy. 
This conservation law was originally derived by Blokhintsev [Q for slowly varying flows, and 
more generally by Cantrell and Hart [Q in their study of the acoustic stability of rocket 
engines. See also reference and |[T^ eq. (5.18). 



Now we turn to the equation for momentum conservation. Working similarly to the 
energy law we find 

= d,T^^ - [Oj, 0]T°° - [zj, 0]r° - [Oj, z]T°^ 

= d^T'^j - p^i)V{^i)idjVi. (5.28) 

Again notice the cancellation of the terms containing p(2). 

The covariant conservation equation D^T'^j = therefore reads 

dtp(^i)V(^i)j + di (^p(o)W(i)i'y(i)j + + 6ijP{2)) + p{i)V{i)idjVi = 0. (5.29) 

The connection terms have provided a source term for the momentum density. Thus, in an 
inhomogeneous flow field, momentum is exchanged beween the waves and the mean flow. 

VI. PHONONS AND CONSERVATION OF WAVE ACTION 



If the mean flow changes only slowly over many wavelengths, the sound field can locally 
be approximated by a plane wave 
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t) = ao cos(k ■ X — a;t). (6.1) 

The frequency uj and the wave-vector k are here related by the Doppler-shifted dispersion 
relation = cj^ + k- v, where the relative frequency, Ur = c\k\, is that measured in the frame 
moving with the fluid. A packet of such waves moves at the group velocity 

]^ 

U = x = c— + v. (6.2) 

\k\ 

As the wave progresses through regions of varying v, the parameters k and oq will slowly 
evolve. The change in k is given by the ray tracing formula (|A16|) 



dt ~ '''dxr ^^-"^^ 

where the time derivative is taken along the ray 

The evolution of the amplitude ag is linked with that of the energy density, Wr through 



{Wr) = ^alp^o)"^. (6.5) 



Now for a homogeneous stationary fluid we would expect our macroscopic plane wave to 
correspond to a quantum coherent state whose energy is given in terms of the (quantum) 
average phonon density N as 

Etot = (Volume) (W^) = (Volume) iV/iu;^. (6.6) 

Since it is a density of "particles", N should remain the same when viewed from any frame, 
consequently the relation 

AT/. = (6.7) 

should hold true generally. In classical fluid mechanics the quantity {Wr)/uJr is called the 
wave action P, p!T|JTB| . 
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The time averages of other components of the energy momentum tensor may be also 
expressed in terms of N. For the mixed tensor we have 

(T°o) = (iy. + v-p(i)V(i)) = Mcu 

(T^o) = ((— + V • V(i))(p(o)t;(i). + p(i)t;.)) = NhuoU, 
P(o) 

(-r°) = (p(i)t;(i).) = Nhh 

{-T)) = {p{Q)V(i)iV(i)j + ViP(i)V(i)j + 5ijP{2)) = NhkjUi. (6.8) 

The last result uses the fact that {P[2)) = for a plane progressive wave. 

Inserting these approximate expressions for the time averages into the Blokhintsev energy 

conservation law ( |5.27| ) we find that 

dNhiu 
~dt 

We can write this as 



+ V ■ {NhuV) = 0. (6.9) 



Nh{^ + V- Vu^j + hio + V • (iVU)j = 0. (6.10) 



The first term is equal to duo/dt along the rays and vanishes for a steady mean flow as a 
consequence of the hamiltonian nature of the ray tracing equations. The second term must 
therefore also vanish. This represents the conservation of phonons, or in classical language, 
the conservation of wave-action. 

In a similar manner the time average of (|5.28|) may be written 



= ^ + V ■ (iVA:,U) + Nk. 



dt \ J J ' I- Q^j 

.^(|l^U.V., + *.|i)..,(f + V.(iVU,). ,6.11) 

We see therefore that the momentum law is equivalent to phonon-number conservation 
combined with the ray tracing equation (|A16|) . 



VII. DISCUSSION 



The possibilty of interpreting the time average of our momentum conservation law in 
terms of quantum quasi-particles should warn us that we are dealing with pseudomomentum 
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and not with newtonian momentum |T8[. Nonetheless the quantity (p(i)V(i)) = NKk is 
rehably computed from the hnearized wave equation, and is part of the true momentum. 
It is simply not all of it. Even in the absence of a mean flow with its (p{2)v) contribution 



we still have to contend with p(o)(v(2)), and this can be important. As an example [jT8 



consider a closed cylinder filled with fluid. At one end of the cylinder a piston is driven 
so as to generate plane sound waves which completely span the cross section of the tube. 
At the other end a second piston is driven at the same frequency with its phase adjusted 
so as to absorb the sound waves without reflection. It easy to see that an extra pressure 
equal to {Wr) is exerted on the ends of the tube over and above whatever isotropic pressure 
acts on the ends and sides equally. It is "obvious" that this is the force per unit area Nhhc 
required to generate and absorb the phonon beam "momentum". Unfortunately for this 
simple idea, it is equally obvious that the time average center-of-mass velocity of the fluid in 
the tube vanishes, so the true momentum density in the beam is exactly zero. The (p(i)V(i)) 
contribution to the momentum density is exactly cancelled by a P(o)(v2) counterflow. This 
eulerian streaming is driven by the fluid source term for (v(2)) implicit in ( p. 161 ) |3CI| ]. (In a 
lagrangian description the particles merely oscillate back and forth with no secular drift). 
The momentum flux however is exactly the same as if (the italics are from |TB[) there was 
no medium and the phonons were particles possessing momentum hh. This is frequently 
true: the flux of pseudomomentum is often equal to the flux of true momentum to O(a^) 
accuracy. Pseudomomentum flux can therefore be used to compute forces. On the other 
hand the density of true momentum in the fluid and the density of pseudomomentum are 
usually unrelated^ 



^This does not mean that the attribution of momentum to a phonon in the two-fluid model for 
a superfluid is incorrect. In superfluid hydrodynamics the /0(o)(v2) counterflow is accounted for 
separately from the (p{i)V(-x)) = iV/ik normal-component mass flux. The counterflow is included 
in the supercurrent needed to enforce V • {pn^n + PsVg) = 0. 
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It should be said that the P(o)(v2) counterflow will not always cancel the P(i)(vi) wave 
pseudomomentum [Q. The (V2) flow depends the geometry. It is found from the source 
equation ( |5.16| ) and from the force the sound field applies to the fluid. The latter will 
be small when there is no dissipation, as is the case in a superfluid, and for an isolated 
sound beam source in an infinite medium (V2) will consist of a flow directed radially inwards 
towards the transducer of sufficient magnitude to supply the mass flowing out along the 



sound beam [^. In the presence of dissipation the force becomes important, leading to 
acoustic streaming. 

Consider our closed cylinder further. From (|4.13|) we see that in a system with fixed (P), 
and in the presence of the sound wave, the mean density of the fluid will be reduced by 

{Wr) rflnc 



(P(2)) 



(7.1) 

P(0) 



d\n.p 

Since our cylinder has fixed volume, this density reduction cannot take place. Instead it is 
opposed by a pressure on the cylinder wall 

dine 



AP = {Wr) 



(7.2) 

P(0) 



d\np 

that must be added to the isotropic pressure in the absence of the sound wave. The complete 
radiation stress tensor is therefore 



This result goes back to Brillouin [^. The true radiation stress therefore differs from the 
pseudomomentum flux tensor in its isotropic part. Forces computed from pseudomomentum 
flux will therefore be incorrect if this pressure gradient is important. Usually it is not. See 



34 1 for examples. 
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APPENDIX A: GEODESICS AND HAMILTONIAN FLOWS 



In this appendix we show that the null geodesies of the Unruh metric coincide with 
conventional hamiltonian optics ray tracing. The usual ray tracing equations are derived 
from a;(k, x) as 

duj ■ duj , _ , 

In our case a;(k, x) = c|A;| + v • k. Thus 

dt \ky dt dx' ^ ^ ' 

We begin by noting that geodesies with an affrne parameter r are stationary paths for 
the lagrangian 

, 1 dx^'dx'' ,^ , 

To make connection with the ray tracing formalism we consider the corresponding hamilto- 
nian 

a = (A4) 

and write down Hamilton's equations with r playing the role of time 

dx^ dH 
dr dpf, 

dp, dH Idg^P 

Combining these gives 
Now for matrices g we have 

dg~' ^ -g-\dg)g-\ (A7) 

so with (g)c,/3 = gap and (g"^)a/3 = g"^ we can write 
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dr^ ^2^ [dx'^ dx^ dx' J dr dr ' ^ ' 

which is the geodesic equation. 

We now examine these equations for the particular case of the Unruh metric. We define 
a 4-vector — (a;, —ki) so that p^x^ = — k • x. Then 

H = Ig^'^P.P. ^Uiu-vkf- c'\k\') . (A9) 



Hamilton's equations become 



dr dr doo 

and 

dx'' dH 



dx^ dt dH , 

— — = — = — = a; - V • k, (AlO) 



Vi{ijj - V • k) + (?ki. (All) 



dr dki 

For null geodesies (a; — v • k)^ — c^j/cp = 0, or (a; — v • k) = c\k\. Thus 



or 



dt \k^ 
which is the group velocity equation. We also find 

duj OH 



dt (a; — V • k) 



^ = + Ct^, (A13) 



dT dt 

if the flow is steady, and 



(A14) 



dT dx^ dx 
which is equivalent to the momentum evolution equation 



dt - dx^"'- ^^^^^ 
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